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ABSTRACT 

Measuring neutron star radii with spectroscopic and timing techniques relies on the combination 
of multiple observables to break the degeneracies between the mass and radius introduced by general 
relativistic effects. Here, we explore a previously used frequentist and a newly proposed Bayesian 
framework to obtain the most likely value and the uncertainty in such a measurement. We find 
that, for the expected range of masses and radii and for realistic measurement errors, the frequentist 
approach suffers from biases that are larger than the accuracy in the radius measurement required 
to distinguish between the different equations of state. In contrast, in the Bayesian framework, the 
inferred uncertainties are larger, but the most likely values do not suffer from such biases. We also 
investigate ways of quantifying the degree of consistency between different spectroscopic measure¬ 
ments from a single source. We show that a careful assessment of the systematic uncertainties in 
the measurements eliminates the need for introducing ad hoc biases, which lead to artificially large 
inferred radii. 

Subject headings: methods: statistical - stars: neutron 


1. INTRODUCTION 

There has been significant recent interest in measuring 
the radii of neutron stars. The radii have been shown 
to be direct probes of the ultradense matter equation of 
state, which, in turn is connected to numerous astrophys- 
ical phenomena, such as the dynamics and outcomes of 
stellar explosions and the signals from the coalescense of 
compact objects. 

Because of strong general relativistic effects, all ob¬ 
servables from the vicinity of a neutron star depend on 
different combinations of its mass and radius. (For fast 
spinning neutron stars, higher order moments such as 
the spin and quadrupole also play a role). As a result, 
measuring radii require at least two distinct observables 
to break the degeneracy with mass. 

During the last decade, measurements of neutron star 
radii have primarily relied on utilizing two spectroscopic 
observables from neutron stars that show thermonuclear 
bursts: the apparent angular size and the Eddington flux 
(Ozel 2006). This technique has been applied to a num¬ 
ber of neutron stars (see., e.g., Ozel et al. 2009; Giiver 
et al. 2010a, b) and led to radius measurements that 
significantly constrained the ultradense matter equation 
of state (Ozel et al. 2010). Similar spectroscopic tech¬ 
niques have also been explored that rely instead on the 
apparent angular size and the Eddington temperature 
(Suleimanov et al. 2012) or the apparent angular size 
and the evolution of the spectral temperature (e.g., Kus- 
mierek et al. 2011). Measurements in the near future 
with NICER (Gendreau et al. 2012) and LOET (Eeroci 
et al. 2012) will utilize different properties of the pulse 
profiles observed from the surfaces of spinning neutron 
stars. In this case, the harmonic content of the pulse pro¬ 
files and their energy dependence provide the necessary 
observables that help break the degeneracy between the 
neutron-star radius and mass (Psaltis et al. 2014). 

In all of these approaches, the measurement of two ob¬ 


servables leads, in principle, to a solution for the two un¬ 
knowns of interest, namely the neutron-star radius and 
mass. In practice, however, the situation is more com¬ 
plicated because of the nonlinear dependence of the ob¬ 
servables on the neutron-star parameters. Depending on 
their particular values, a set of measurements may lead 
to two, one, or no solutions for the mass and radius. 
Because of these nonlinearities, the inferences from such 
measurements likely depend on the particular statistical 
estimators used in combining these constraints. 

In this paper, we assess, using mock data, the previ¬ 
ously used frequentist method for measuring neutron star 
radii based on two spectroscopic observables. We show 
that the radii inferred with this frequentist approach are 
often biased and, in some cases, the true solution is for¬ 
mally excluded. The level of bias increases significantly 
as the measurement errors increase. We then devise a ap¬ 
proach within a Bayesian framework that is substantially 
less affected by these nonlinearities, even in the case of 
large uncertainties. We also discuss the conditions un¬ 
der which the lack of solutions for a large fraction of 
the parameter space can be used to infer the presence of 
systematic uncertainties in the measurements. 

2. THE PREVIOUS APPROACH TO DETERMINING RADII 

We will focus hereafter on the measurement of the neu¬ 
tron star radii and masses based on the combination of 
the apparent angular size and the Eddington flux ob¬ 
served during thermonuclear bursts. Even though we 
use this as our primary example, the results are general 
and can be translated to the other combinations of ob¬ 
servables discussed in the introduction. 

We follow the formalism in Ozel et al. (2009) and as¬ 
sume that the apparent angular size A and the Edding¬ 
ton flux Ftd have been measured during the cooling tails 
and the touchdown moments of thermonuclear bursts, 
respectively. In the Schwarzschild approximation, which 
is appropriate for slowly spinning neutron stars, these 
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quantities are related to the neutron star mass M and 
radius R according to 


and 
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where G is the gravitational constant, c is the speed of 
light, D is the distance to the neutron star, fc is the color 
correction factor due to the stellar atmosphere, fees = 
0.2 (1+X) cm^ g“^ is the electron scattering opacity, and 
X is the hydrogen mass fraction. Note that for moder¬ 
ately spinning stars, there are additional corrections that 
depend on the spin and the quadrupole moment of the 
neutron star (see Baubock et al. 2015). In addition, the 
Eddington flux is subject to temperature-dependent cor¬ 
rections due to the energy-dependent terms in the Klein- 
Nishina cross section, as discussed in Paczynski (1983). 
For simplicity, we ignore these corrections in the present 
statistical treatment; see, however, Ozel et al. (2015) for 
their effects on the inferred radii and masses. 

In principle, these two equations can be solved for the 
neutron star mass and radius given the observables. Be¬ 
cause of the nonlinear nature of these equations, there 
can be zero, one, or two solutions. The number of so¬ 
lutions depends on the value of the quantity (see, e.g., 
Steiner et al. 2010) 


Ftd^esR 
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such that, when a > 1/8, there are no solutions, when 
a = 1/8, there is one (double) solution, and when a < 
1/8, there are two distinct solutions. The critical value 
occurs when R = AGM/c^ independent of all the other 
parameters. 

When this method is applied in practice, the uncertain¬ 
ties inherent in the measurements of A and Ftd need to be 
converted into uncertainties in R and M . In this frequen- 
tist approach, this is achieved by sampling the likelihoods 
over the observables, calculating the mass and radius for 
each pair of observables, and using this to populate a 
posterior likelihood over radius and mass. In Ozel et al. 
(2009), this was carried out analytically using the Jaco¬ 
bian transformations of the posterior likelihoods. Steiner 
et al. (2010) repeated the same frequentist analysis using 
Monte Carlo techniques to sample the likelihoods over 
the observables. Both methods of calculation give iden¬ 
tical results under the same set of assumptions. 

In this formalism, the posterior likelihood over mass 
and radius is given by 

P{M,R)dMdR=^ j P{D) dD J P{f,) dfcj P{X) dX 
xP[Ftd{M, R, D)]P[A{M, R, D)] J dMdR. (4) 


Here, P{D), P{Ftd), and P{A) are the likelihoods over 
the distance, touchdown flux, and the apparent angular 
size measurements, respectively. P{X) and P{fc) are 
the priors over the hydrogen mass fraction and the color 


correction factor. Finally, 
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is the Jacobian of the transformation. The factor 1/2 
appears in this equation because nearly all pairs of ob¬ 
servables correspond to two distinct pairs of M and R. 
The region of the parameter space for which the observ¬ 
ables correspond to one repeating solution for M and R 
has zero volume and will, therefore, not contribute to the 
final likelihood. 

It is evident from equations o and dD) that P(M, R) is 
identically equal to zero when R = AGM/c^, independent 
of the measurements. This occurs because for masses 
and radii that satisfy this condition, the two constraints 
imposed by the observables (eqs. [T] and [5]) are parallel 
at the point of contact and the system is degenerate. 
As a result, if the mass and radius of the neutron star 
he along the R = AGMjt? line, which is very likely for 
the expected range of neutron star masses and radii, the 
posterior likelihood will exclude the true solution and, 
hence, will introduce a bias in this case. 

To explore the severity of this bias for a range of neu¬ 
tron star masses and radii, we perform the frequentist in¬ 
ference described above for three different pairs of mock 
measurements that correspond to stars with different 
RjM ratios. In all cases, we calculate the apparent an¬ 
gular size and the Eddington flux for the assumed mass 
and radius of the neutron star by fixing the hydrogen 
mass fraction to A = 0, the distance to 4 kpc, and the 
color correction factor to 1.35. We then assign a 5% 
Gaussian error in A and Ptd- To infer R and M from 
these mock measurements, we assume a perfect prior 
knowledge of the hydrogen mass fraction and distance 
but take a boxcar prior in the color correction factor be¬ 
tween 1.3 and 1.4. The left panels of Figure [T] show the 
68% and 95% confidence contours of the posterior like¬ 
lihood over the inferred mass and radius while the right 
panels show the posterior likelihood over radius when 
P(i?,M) is marginalized over mass. In the left panels, 
the black lines correspond to contours of constant A and 
Ftd at their central values and for /c = 1.35, such that 
one of the intersection of these two lines correspond to 
the assumed radius. The green line is the critical curve 
R = AGM/c^. 

The examples depicted in Figure 1 show that when 
the assumed value of i?/M is significantly away from the 
critical curve, the posterior likelihood of this frequen¬ 
tist approach is centered on the true solution, with little 
evidence for bias. However, as the assumed ratio ap¬ 
proaches the critical curve, which also corresponds to 
more realistic values of masses and radii, the inferred 
radii can be substantially biased toward higher or lower 
values, by as much as 1.5 km. This is an unacceptably 
large bias given that ^ 1 km precision is required in the 
radius measurements in order to place meaningful con¬ 
straints on the equation of state (see, e.g., Ozel & Psaltis 
2009). Moreover, in this idealized example, we assumed 
that the distance is known with extremely high accuracy, 
which is unrealistic for the neutron star sources used for 
the radius measurements. In Section 4, we will explore 
the effect of realistic distance uncertainties on the radius 
measurement bias and show that it gets even worse in 
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this frequentist approach. For the measurenients per¬ 
formed to date using this approach, this implies that the 
inferred radii are likely to have been overestimated by 
-1-1.5 km. 

3. A NEW BAYESIAN APPROACH TO DETERMINING 

RADII 

The posterior likelihoods over the neutron star mass 
and radius can be inferred from the same observables 
P{A) and P{F^^) in a different way using Bayes’ theorem, 
according to which 

P(M, i?|data) = C'P(data|M, R)Ppri(M)Ppri(R), (6) 

where Pp^i{M) and Ppri(i?) are the priors over the mass 
and radius and C is an appropriate normalization con¬ 
stant. Given that A and Ftd are ideally uncorrelated 
measurements, we can write 

P(data|M,R) = J P{D) dD J P{f^) dfc J P{X) dX 
xP[FtdiM,R,D)]P[AiM,R,D)]. (7) 

Assuming a flat prior over the radius and mass, it is clear 
from equations (|3]) and © that the difference between 
the Bayesian and frequentist inferences discussed here is 
the presence of the Jacobian in the latter. Consequently, 
the biases introduced in the frequentist approach by the 
zero points of the Jacobian are not present here. The 
underlying reason is that in a Bayesian framework, one 
does not ask what is the mass and radius pair that corre¬ 
sponds to each set of observables, but rather, for a given 
mass and radius, what is the chance of obtaining the 
corresponding set of observables. 

In the left panels of Figured! we overplot the 68% and 
95% confidence contours of the posterior likelihood cal¬ 
culated according to equation [7| and show the marginal¬ 
ized likelihood over the radius in the right panels. When 
the R/M ratio is significantly away from the critical 
curve, the Bayesian approach gives unbiased results as 
in the frequentist approach. There is still a difference 
between the two approaches, however, which becomes 
evident in the marginalized likelihoods: the Bayesian ap¬ 
proach gives equal maximum likelihood between the two 
equivalent solutions, whereas the frequentist yields equal 
integrated likelihoods in the two islands. Given that the 
two solutions are mathematically equivalent (albeit the 
smaller radius is sometimes physically unacceptable), the 
difference in the integrated posterior likelihoods between 
the two islands should not be used to discriminate be¬ 
tween the two solutions. 

The biggest difference in the two frameworks occurs 
when the R/M ratio of neutron stars approach or lie 
on the critical curve. In these cases, the 68% confi¬ 
dence contour of the Bayesian posterior likelihood en¬ 
compasses both solutions and, both in the full P{M, R) 
and in the marginalized P{R) likelihoods, the method re¬ 
sults in larger uncertainties in the inference of the radii 
than in the frequentist approach. However, it does not 
suffer from the biases of the frequentist approach, which 
are exacerbated when the neutron stars lie on the critical 
curve (see the bottom panels of Fig. [IJ. 

4. BIASES IN THE INFERRED EQUATION OF STATE 

FROM MULTIPLE RADIUS MEASUREMENTS 


In the previous sections, we addressed the biases in the 
inferred radius of a single star using two approaches, a 
previously developed frequentist and a newly proposed 
Bayesian formalism, under the assumption of small un¬ 
certainties in the two measured quantities A and Ftd. 
In reality, multiple radius measurements, with different 
sources of uncertainty, are necessary to place significant 
constraints on the dense matter equation of state. In this 
section, we compare the performances of these two ap¬ 
proaches under frequentist considerations. In particular, 
we explore the uncertainties and biases in the inferred 
properties of the equation of state using mock data sets 
for multiple objects. 

Earlier work has shown that all equations of state that 
smoothly connect to the low density equation of state 
of matter can be represented by a small number of pa¬ 
rameters, which can be inferred by the measurement of 
several neutron star radii at different masses (Lindblom 
1992; Read et al. 2009; Ozel & Psaltis 2009; Gandolfi 
et al. 2014). In the relevant range of observed neutron- 
star masses (see, e.g., Ozel et al. 2012), most equations 
of state predict nearly constant radii. For simplicity in 
the present study, we make use of this fact to reduce the 
complexity of the equation of state and represent them 
by a single parameter, which we take to be the constant 
radius. 

For each assumed neutron star equation of state, and, 
hence, assumed radius, we draw five neutron star masses 
from a Gaussian distribution with a mean of 1.46 Mq 
and a dispersion of 0.21 Mq; this is the distribution of 
masses that is observationally inferred for the descen¬ 
dants of X-ray bursters (Ozel et al. 2012). For each of 
these neutron stars, we calculate the Eddington flux and 
the apparent angular size assuming a distance of 4 kpc, 
color correction factor of 1.35, and a hydrogen mass frac¬ 
tion X = 0. We assign a 10% Gaussian uncertainty in 
each of these values, which is typical for X-ray bursters 
(see Giiver et al. 2012 a,b) and draw randomly a set of 
mock measurements from these distributions. We then 
apply the frequentist and Bayesian frameworks discussed 
above to each of these measurements, assuming that the 
distance and the hydrogen mass fraction are known a 
priori, but allow the color correction factor to be in the 
range 1.3-1.4. We marginalize the five resulting posterior 
likelihoods over mass and then multiply them so that we 
obtain the posterior likelihood over the inferred radius 
for that equation of state. We repeat this procedure for 
10,000 equations of state with assumed radii between 9 
and 14 km and different realizations of the neutron star 
masses drawn from the same distribution. 

In Figured! we plot the inferred vs. the assumed ra¬ 
dius for each set of mock data for the frequentist and the 
Bayesian inferences. As expected from the discussion 
in Section 2, the previously used frequentist inference 
yields two well separated tracks of solutions, whereas the 
proposed Bayesian approach has a single band. (Note 
that the few outliers correspond to the unlikely situa¬ 
tions where all five drawn masses are comparable and 
place the neutron stars away from the critical curve). 
More importantly, even if the band of small radii is re¬ 
jected on physical grounds, the remaining band in the 
frequentist approach yields measured radii that are bi¬ 
ased toward larger values. In contrast, the radii inferred 
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Fig. 1.— (Left) Contours of 68% and 95% posterior likelihoods over mass and radius using (red) the frequentist and (blue) Bayesian 
approaches discussed in the text, for three neutron stars with different assumed masses and radii. The true masses and radii are denoted 
by red filled circles. The black curves correspond to lines of constant apparent angular size, A, and Eddington flux, The green line 

corresponds to the critical curve R = AGMfc^. (Right) The same posterior likelihoods marginalized over mass. As the true mass and radius 
of a neutron star approaches the critical curve, the frequentist approach significantly biases the inferred radii. In all cases, the distance of 
the neutron star was fixed to 4 kpc, the color correction factor to fc = 1.35, and the hydrogen mass fraction to A = 0. 
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Fig. 2.— The radius inferred for an equation of state (chosen 
here to predict constant radius neutron stars over the mass range 
of interest) using mock spectroscopic observations of five neutron 
stars, plotted against the assumed radius that it predicts, for a 
large ensemble of realizations of the mock data, in the (upper) fre- 
quentist and (lower) Bayesian approaches discussed in the text. 
All simulated neutron stars were placed at a distance of 4 kpc 
and their masses were drawn from a Gaussian distribution that 
has been inferred observationally for the decendendents of X-ray 
bursters. The spectroscopic measurements were assumed to have 
10% Gaussian errors, while the distance and hydrogen mass frac¬ 
tion was assumed to be accurately known a priori. The frequentist 
approach leads to substantial biases in the infered radii, while the 
Bayesian approach does not suffer from biases. 

in the Bayesian formalism are centered on the assumed 
values, with little, if any, hints of bias. 

To quantify the degree of bias in these measurements, 
we plot in Figure [3] the histogram of the difference be¬ 
tween the inferred and assumed radii for two different 
values of the assumed radius. The bias in the frequen¬ 
tist approach (again focusing on the larger of the two 
solutions) is -1-0.75 km for a 12 km neutron star and 
-hi km for a 10 km neutron star. The histogram in the 
Bayesian approach, on the other hand, is centered on the 
true value for both assumed radii, even though it is wider 
than either of the frequentist peaks. In other words, the 
frequentist approach generates the more precise but less 
accurate solution. 

The last aspect of the measurements that may bias 
the inferred radii is how well the distance to the source 



Radius Measured - Radius Assumed (km) 



Radius Measured - Radius Assumed (km) 

Fig. 3.— The histogram of the difference between the inferred 
and assumed radii for two different values of the assumed radius, 
for the simulated data shown in Figure |2] The bias toward larger 
radii in the frequentist approach is ~ 0.75 km for 12 km neutron 
stars and ~ 1 km for 10 km stars. The proposed Bayesian inference 
is more accurate but less precise. 

is known, as reflected in the width of the prior over the 
distance P{D) that is considered in the calculation. This 
is of particular concern in the frequentist approach be¬ 
cause the Jacobian of the transformation scales as D* 
and, therefore, places large weight on the smaller dis¬ 
tances. We repeated the above simulation of mock ob¬ 
servations with two differences. First, we reduced the 
assumed uncertainties in the measurements of A and Ftd 
to 5% in order to focus on the effect of the broader range 
of distances. Second, even though we placed all simu¬ 
lated neutron stars at 4 kpc, we assumed a flat prior in 
the distance in the range 3-5 kpc when inferring the radii 
from the mock measurements. 

In Figure 01 we plot the inferred vs. the assumed ra¬ 
dius for each set of mock data when the range of distances 
in the priors dominate the uncertainties of the measure¬ 
ment. As in the previous set of mock data, there are 
two bands of acceptable solutions. However, we will not 
consider the bands that correspond to small radii any 
further, as these can typically be rejected on physical 
grounds. In the frequentist inference, the large range 
of distance priors biases the inferred radii toward values 
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Fig. 5.— (Left) The expected distribution over the parameter a 
calculated for the observed distribution of neutron masses and for 
several values of the neutron star radius that are within the physi¬ 
cally reasonable range. The possible values of a. span an extremely 
narrow range and the likelihood peaks sharply toward the critical 
value a. = 1/8 for 10 — 12 km neutron stars. 


tance D and the hydrogen mass fraction X are flat over 
a large range of values. As a result, the range of inferred 
a values span a wide range, with only a fraction of them 
falling below the 1/8 limit. This fraction was used by 
Steiner et al. (2010) to argue for the inconsistency of 
observables from some sources; the same argument was 
later repeated by Miller (2013). We will now use mock 
data to explore this issue quantitatively and will show 
the shortcomings of this argument. 

We first calculate the expected distribution of a values 
for realistic neutron star masses and radii. Using equa¬ 
tions (1), (2), and (3), we can write a only in terms of 
M and R as 


Fig. 4.— Same as Figure [2] but for assumed 5% uncertainties 
in the spectroscopic measurements and a flat prior in distances in 
the range 3-5 kpc. In this case, the frequentist approach shows a 
significant bias toward lower radii, while the Bayesian approach is 
only marginally affected. 

that are smaller by as much as 1 km than the assumed 
radii. The reason for this bias is the strong dependence 
of the Jacobian on distance, which favors the smaller 
distances and hence the smaller radii. In contrast, the 
Bayesian inference suffers from a much smaller bias to¬ 
ward smaller radii, even in this case. 


GM / 2GM\ 
~ F? V ~ i?c2 ) ■ 


( 8 ) 


If we assume that all neutron stars have the same radius 
i?o in the mass range of interest and masses drawn from 
the observed distribution of their descendants (i.e., fast 
pulsars, see §2), then we can write the expected distri¬ 
bution over a as 


P{a) = P{M) 


dM 

da 


(9) 


5. THE POSTERIOR PROBABILITIES OF THE 
CONSISTENT MASS-RADIUS SOLUTIONS 

In the previous sections, we investigated the biases and 
uncertainties in the spectroscopically inferred masses and 
radii of neutron stars within a frequentist and a Bayesian 
framework. We showed that the fact that real solutions 
for mass and radius can be obtained for sets of observ¬ 
ables that satify a < 1/8 (see eq.|3|) introduces significant 
biases in this frequentist approach, especially when the 
true M, R of the neutron star lies near this critical value. 

We now turn our attention to a second concern aris¬ 
ing from the presence of this limit on a. In principle, 
the situation where a set of observables with negligible 
uncertainties yield a > 1/8 would imply the inconsis¬ 
tency of the observables and priors with one another. In 
practice, however, the two observables Ftd and A have 
uncertainties and, more importantly, the priors over dis- 


Equation ([5]) has two solutions 


_ 

M± = — {i±vr^) ( 10 ) 

and, therefore, for a < 1/8, 


2(t^ I 2(T^ 

e M + e M 

( 11 ) 

Here, Mq = 1.46 Mq and um = 0.21 Mq (Ozel et al. 
2012). We plot in Figure [5] the distribution over a for 
several physically reasonable neutron star radii. 

For radii in the 10-11 km range, the expected distribu¬ 
tion over a is very narrowly peaked, with an integrable 
pole at the limiting value of a = 1/8. Asa result, for such 
equations of state, obtaining a values from observations 


P(a) = 


G\J\ — 8q! y27rcr 
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that are close to the critical value is, in fact, expected, 
and cannot be taken as an indication for the inconsis¬ 
tency of the observables, as suggested by Steiner et al. 
(2010). More importantly, because the two observables 
Ftd and A have non-negligible uncertainties and the pri¬ 
ors over D and X are typically flat over a wide range 
of values, the observationally inferred distribution of a 
will be substantially broader and will very often peak at 
values a > 1/8. It then follows that a very large section 
of the observationally inferred distribution over a will be 
rejected as leading to unphysical solutions even for fa¬ 
vorable uncertainties, without implying an inconsistency 
among the observables and the priors. 

In order to demonstrate this point, we simulate mock 
data for a 1.7 Mq, 10 km neutron star placed at a dis¬ 
tance of 5 kpc, with a hydrogen mass fraction X = 0 dur¬ 
ing Eddington-limited bursts and an atmospheric color 
correction factor /c = 1.35. We assign Gaussian 2% un¬ 
certainties in the apparent angular size A and in the 
Eddington flux Fta; these are comparable to the for¬ 
mal uncertainties reported by Giiver et al. (2010a) for 
4U 1608—52. We then draw Monte Carlo pairs from 
these two distributions as our mock data for these two 
quantities. We also take a mock measured distance that 
is flat between Umin and Umax such that the true dis¬ 
tance lies between these two values. We also assume, as 
is sometimes the case, that there is no measurement of 
the hydrogen mass fraction and we treat is as a flat prior 
between W = 0 and X = 0.7. 

The fraction of “accepted” solutions that correspond 
to q; < 1/8 can be written as an integral over the entire 
parameter space, limited by the condition 




3y2^1/2 


Sko{l+X)D 


( 12 ) 


such that 
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.(13) 


Here, Hobs and E/bs are the mock observed data for each 
realization and we set the observed uncertainties a a and 
ap equal to the assumed ones; i.e., the measurements 
do not over- or underestimate the true uncertainties. In 
the top panel of Figure [SI we show the fraction ^ of ac¬ 
ceptable solutions as a function of the upper limit that 
can be placed on the source distance; for these simula¬ 
tions, we set Emin = 4.9 kpc. Even in the unrealistic 
case where the distance is extremely well constrained, 
i.e., Emax — Emin < 0.2 kpc, the fraction of accepted 
solutions is equal to at most a few percent. As the up¬ 
per bound on the distance increases, the fraction drops 
to 0.1% or less. Such values are comparable to the ones 
obtained by Steiner et al. (2010) for 1608—52, on which 
these simulations were based, and show that fractions of 
this order clearly do not imply any level of inconsistency. 

In the bottom panel of Figure |6l we show the results 
of a similar simulation in which the observational uncer¬ 
tainties are underestimated by a factor of 5; i.e., the dis¬ 
persions of the distributions used to draw the mock data 
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Fig. 6.— The fraction of consistent solutions in the frequentist 
approach as a function of the maximum allowed distance in the 
distance prior. (Top) The mock data are generated for a 1.7 Mq, 
10 km neutron star assumed to be at 5 kpc and with hydrogen mass 
fraction X = 0. They are drawn from Gaussian distributions with 
2% uncertainty in in Ftd and A, centered at the true values for the 
assumed parameters. In the calculation of the fraction of consistent 
solutions, the priors in the distance and the hydrogen mass fraction 
are assumed to be flat, spanning the range from 4.9 kpc to Umax 
and 0.0 to 0.7, respectively. The fraction of consistent solutions 
(see the text for definition) is typically ~ 1% is known accurately 
a priori and is dramatically reduced to < 0.1% as the upper bound 
on the source distance increases. (Bottom) Same as above but with 
true uncertainties in the observables taken to be five times larger 
than the assumed uncertainties. Underestimating the uncertainties 
by a factor of 5 reduces the fraction of consistent solutions to as 
low as 10~®. 

are 10%, while a a and ap are still set to 2%. We are 
exploring this situation because the subsequent compre¬ 
hensive analysis of all the X-ray data from sources that 
show a large number of thermonuclear bursts indicates 
a ~ 10% spread in the inferred values of A and Ftd for 
each source, which is larger than the formal uncertain¬ 
ties of the individual measurements (Giiver et al. 2012a, 
b). In this case, the fraction drops substantially. Steiner 
et al. (2010) used the similarly low values inferred for 
4U 1820—30 when only the statistical uncertainties in the 
data were considered to justify an ad hoc and unphysical 
reinterpretation of the Eddington limit. This increased 
the fraction of consistent solutions by artifically moving 










neutron stars to larger radii and thus, away from the 
critical value of a = 1/8. Instead, a correct assessment 
of the systematic uncertainties in the spectroscopic mea¬ 
surements, as was done in the later studies, provides a 
much simpler explanation of the inferred small fraction 
of accepted solutions, without the need for ad hoc rein¬ 
terpretations of the data. 

6. CONCLUSIONS 

Measuring neutron star radii with spectroscopic and 
timing techniques requires at least two observed quanti¬ 
ties to break the degeneracies between the mass and ra¬ 
dius introduced by general relativistic effects. The sys¬ 
tem of equations that connect the observables to mass 
and radius often have critical points and regions of no 
solutions. In this paper, we explored the biases in the 
inferred radii introduced by this mathematical property 
of the problem. 

We assessed a previously used frequentist method for 
inferring neutron star radii, devised a new one within a 
Bayesian framework, and compared their performances 
under frequentist considerations. We showed that the 
former suffers from significant biases in the range of 
masses and radii predicted by the modern equations of 
state and when realistic uncertainties are taken into ac¬ 
count in the measurement of the distance to the source. 
In contrast, in the latter framework, the inferred uncer¬ 
tainties are larger but the most likely values do not suffer 
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from such biases, making it the preferred framework for 
measuring neutron star radii using current and future 
data. 

Finally, we explored ways of quantifying the degree 
of consistency between different spectroscopic measure¬ 
ments from a single source. We demonstrated that the 
fraction of the parameter space over the observables that 
gives rise to real solutions in mass and radius is not 
a good measure of the consistency, as was previously 
claimed. Fractions of accepted solutions of order 0.1% 
are common due to the intrinsic uncertainties in the dis¬ 
tance and hydrogen mass fraction. In addition, signif¬ 
icantly smaller values of the accepted fraction can be 
accounted for by the recently determined systematic un¬ 
certainties in the spectroscopic measurements. In a com¬ 
panion paper, we will apply the statistical framework 
discussed here to all the currently available spectroscopic 
measurements of the neutron star radii and constrain the 
dense matter equation of state. 
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